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ABSTRACT 

We have performed VLBI and single-dish monitoring of 22 GHz H2O maser emission 
from the high mass young stellar object G353. 273+0.641 with VERA (VLBI Explo- 
ration of Radio Astrometry) and Tomakamai 11-m radio telescope. Two maser flares 
have been detected, separated almost two years. Frequent VLBI monitoring has re- 
vealed that these flare activities have been accompanied by structural change of the 
prominent shock front traced by H2O maser alignments. We have detected only blue- 
shifted emissions and all maser features have been distributed within very small area 
of 200 x 200 au 2 in spite of wide velocity range (> 100 km s" 1 ). The light curve shows 
notably intermittent variation and suggests that the H2O masers in G353. 273+0. 641 
are excited by episodic radio jet. The time-scale of ~ 2 yr and characteristic velocity 
of ~ 500 km s _1 also support this interpretation. Two isolated velocity components of 
C50 (-53 ± 7 km s" 1 ) and C70 (-73 ± 7 km s" 1 ) have shown synchronised linear ac- 
celeration of the flux weighted Vlsr values (~ -5 km s _1 yr -1 ) during the flare phase. 
This can be converted to the lower-limit momentum rate of 1.1 x 10~ 3 Mq km s" 1 
yr" 1 . Maser properties are quite similar to that of IRAS 20126+4104 especially. This 
corroborates the previous suggestion that G353. 273+0. 641 is a candidate of high mass 
protostellar object. The possible pole-on geometry of disc-jet system can be suitable 
for direct imaging of the accretion disc in this case. 

Key words: ISM: jets and outflows - masers - stars: early-type - stars: formation. 



1 INTRODUCTION 

Astronomical masers are variable phenomenon seen in a high 
mass star-forming region. Interstellar H2O masers (Jk lx k c = 
616-523) at 22. 23508GHz, which appear at the earliest phase 
of formation, are highly variable in particular. They are of- 
ten excited in strong shocks on a working su rface between 
proto s tellar outflow and dense envelope (e .g., iFuruva et al.l 
l2000l : iMoscadelli, Cesaroni fc Rioial l2000l h Outflow prop- 
erties such as size, velocity and morphology, predicted 
from interferometric observations of H2 O masers vary 
source to source (e.g., collimated jet (Shepherd et al. 
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l2004lh wide angle flow (|Motogi et al.ll2008l'). expanding shel l 
l|Torrelles et~atll2003h . e quatorial flow ((Trinidad et alj|2007h 



or combination of them |Torrelles et all 20 111 )'). It is possibly 
related to a number of factors such as stellar mass, evolu- 
tional stage, geometry of surrounding envelope and driving 
mechanism of host outflow. Several single-dish based stud- 
ies h ave reported su b-yr scale maser flares in some cases 
fe.g.. lFelli etai1l2007r i. This time-sale is comparable with or 
even shorter than that of episo dic radio jet, w hich traces a 
root of protostellar outflow (e.g.,[Xnglada 199(| and reference 
therein), and indicates that some portion of H2O masers 
contains valuable information about an innermost outflow. 
Exact relation between an outflow activity and maser vari- 
ability is still unclear, but once it is set up, H2O maser will 
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be an excellent tool to survey dynamic outflow activities in 
small scale, since its bright emission allows us easy and fre- 
quent monitoring even with a small size radio telescope. 

In this paper, we report on the new VLBI and single- 
dish study for intermittent flare activities of H2O masers 
in G353. 273+0.641 (hereafter G353). G353 is a strong 22 
GHz H2O maser s ite in the southern h igh mass star- forming 
region NGC6357 dSakellis et al.lll984l) . The source distance 
is 1.7 kpc from the sun l|Neckellll978l). Multi-epoch ATCA 
observation has been reported in lCaswell k, Phillips! (|2008T ) 
(hereafter CP08). Class II CH 3 OH maser emission [Jk — 5i— 
60 A + ) at 6.668519 GHz is also associated with this source 
(CP08), and hence, the host young s tellar object (YSO ) 
is identified as a high mass YSO (e.g., [M 

uiicr et al.ll2003h . 

CP08 has also suggested that G353 is still in the pre-ultra 
compact (UC) HII region phase, i.e., high mass protostellar 
phas e, based on the absence of any detectable OH masers 
fe.g.. lCaswellll997l ; lBre~en et alj|20ld ). 

G353 has been classified as dominant blue-shifted H2O 
maser in CP08. That is, almost all flux is concentrated on 
blue-shifted emission, despite of very broad velocity range 
of ± 100 km s" 1 with respect to the systemic velocity of 
-5 km s _1 . They have argued that this type of masers can 
be caused by well-collimated jet aligned close to the line of 
sight. There are a few maser sources which show similar blue- 
shift dominance re.g.. lCasweill2004l : CP08: ICaswell fc Breenl 
2010). Some of them suggest acceleration of outflowing ma- 
terials in the jet. The statistical analysis in lCaswell fc Breenl 
(2010) indicates that such a blue-shift dominance is a char- 
acteristic of H2O masers at the earliest evolutionary stage 
of star-formation. 

Although the signature of yr-scale variability of the H2O 
maser has been pointed out by CP08, their small dataset is 
still insufficient for detailed discussion. Our VLBI and single- 
dish monitoring with VERA (VLBI Exploration of Radio 
Astrometry) and the Hokkaido University Tomakomai 11- 
m radio telescope in Japan, which has been started from 
November 2008, has provided more tight-interval dataset 
and detected intermittent flare activities accompanied by 
spatial changes of the shock front traced by maser distribu- 
tion. 



were divided into 512 spectral channels yielding final spec- 
tral resolution of 0.21 km s _1 . 

Special calibration methods for VERA 's dual-beam 
instru ment have already be en described in Motogi et al.l 
|201ll ). We used J1717-3342 jPetrov et al.ll2006l ) as a phase 
calibrator. This source is separated from G353 by 1°.83 at 
a position angle of -108° east of north. J1717-3342 has been 
enough bright (~ 0.7 Jy) and shown no significant struc- 
ture during two years . We also observed NRAO530 (=J1733- 
1304: iMa et akll 19981 ) as a delay and bandpass calibrator ev- 
ery 120 minutes. We have not observed any flux calibrator 
and antenna gains have been corrected based on measured 
system noise temperatures (T sys ) of each station. Typical 
T sys values were ~ 350 K for Mizusawa and Iriki stations 
and ~ 500 K for Ogasawara and Ishigaki stations in the case 
of normal weather condition. Total on-source time is about 
2.5 hour in the first five epochs and 5 hour in the last three 
epochs. We note that the observations were time-shared with 
another maser - calib rator pair, G5. 89-0. 39 - J1755-2232, in 
the first five epochs (|Motogi et al.ll201l[ ). This causes rela- 
tively short on-source time in the former five epochs, but 
total hour angle range (i.e., achieved UV - coverage) is al- 
most same as the latter three epochs. 

Data reduction was carried out using NRAO Astro- 
nomical Imaging Processing System (AIPS) package. We 
adopted standard r eduction procedure for VERA data (See 
iMotogi et al.l l201ll in details) , where we have first deter- 
mined absolute coordinates of several maser spots in 2- 
beam phase referencing analysis, and then, performed self- 
calibration using the brightest maser spot in each epoch. 
We finally searched all maser spots within l".5xl".5 area 
with 7-<t detection limit. Here, the term maser spot indi- 
cates a maser component seen in a single velocity channel, 
and hereafter, the term maser feature means physical gas 
clump which is consist of several maser spots detected in 
successive velocity channels and closely located each other. 

Table 1 summarised all the VERA observations. It 
contains observing dates, synthesised beam sizes in milli- 
arcsecond (mas), beam position angles (PA) east of north, 
typical noise levels and dynamic ranges of synthesised im- 
ages for the brightest maser spot, estimated errors of abso- 
lute positions (see section 3.1) and brief comments, if any, 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 VLBI observations 

VLBI monitoring o f H2O masers with VERA 
|Kobavashi et all 120081 ) has been performed eight epochs 
from November 2008 to April 2010 and is still on-going. 
Each observation was made in VERA's dual-beam mode 
in which targeted maser sou rce and phase calibrator were 
observed simultaneous ly (Ka waguchi. Sasao fc Manabd 
l2000l : lHonma et al.ll2003h . 

Left-handed circular polarised signals were quantised 
at 2-bit sampling and filte red with the VERA digital fil- 
ter unit ijlguchi et al.ll2005h . after that, data were recorded 
onto magnetic tapes at a data rate of 1024 Mbps. The data 
correlation was perf ormed with the Mitaka FX correlator 
ijShibata et al.l 1 19981 ). Total bandwidths assigned for cali- 
brators and the maser are 240 MHz (16 IFs) and 8 MHz 
(1 IF), respectively. Whole velocity coverage of 108 km s -1 



2.2 Single-dish observations 

Single-dish monitoring has been done using th e Hokkaido 
Univ ersity Tomakomai 11-m radio telescope (|Sorai et al.l 
2008). Observing interval is about 1-3 times per month 
except for the receiver maintenance session during July to 
October 2009. The beam width (FWHM) and aperture effi- 
ciency is 4'. 2 and 49 % at 22 GHz, respectively. These values 
h ave been slightly im proved compared with that mentioned 
in lSorai et al.l l|2008l ) because of the upgrade in optics. The 
pointing accuracy is better than 30" in each of azimuth and 
elevation direction. Observations were performed in the po- 
sition switching method. We received single left circular po- 
larisation with 16-MHz bandwidth (216 km s~ in velocity) 
and it was divided into 2048 spectral channels. An intensity 
scale was calibrated by a chopper wheel method every 10 
minutes. 

Observed data were reduced with the NEWSTAR soft- 
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Table 1. Summary of VERA Obscvations 



Epoch 


Date 


Synthesised Beam 


PA 


1-ct Noise" 


DN p b / DN S 


Positional Error 


Comments 






(mas X mas) 


(°) 


(Jy beam -1 ) 




(mas X mas) 


for T sys 


1 


2008/7/1 


3.16 X 0.74 


-17.1 


0.48 


13 / 117 


0.15 X 0.49 




2 


2008/11/18 


3.38 x 0.78 


-21.7 


0.30 


6 /49 


0.92 x 2.15 


~ 1000 K at OG 


3 


2009/2/6 


2.63 x 0.91 


-14.1 


0.37 


14 / 351 


0.06 x 0.24 




1 


2009/5/18 


2.57 x 0.88 


-13.2 


0.30 


16 / 1122 


0.04 x 0.16 




5 


2009/9/16 


2.86 x 0.82 


-12.3 


0.37 


10 / 19 


0.11 x 0.48 




6 


2009/11/5 


2.79 x 0.79 


-13.3 


0.16 


15 / 338 


0.23 x 0.98 




7 


2010/1/21 


3.15 x 0.85 


-19.5 


0.21 


15 / 72 


0.09 x 0.25 


~ 2000 K at OG 


8 


2010/4/21 


2.59 x 0.75 


-14.1 


0.40 


12 / 184 


0.10 x 0.41 


~ 1500 K at IS 



a Typical value in self-calibrated images. 

b The dynamic ranges of the brightest spot in phase- referenced (DN P ) and self-calibrated image (DN S ). 
c IR, IS, OG:Iriki, Ishigaki, Ogasawara station, respectively. 



ware package developed in Nobeyama Radio Observatory. 
Initial 2048 spectral channels were binned up every two 
channels, in order to obtain the same spectral resolution 
as VERA data. Typical rms noise level is ~ 3 Jy in winter 
and increases up to 6 - 25 Jy in the summer season, 



3 RESULTS 

3.1 Spatial and velocity distribution 

Figure 1 shows the spatial and velocity distribution of H2O 
maser spots in G353. All detected maser spots in eight 
epochs are exhibited in this figure, since several maser fea- 
tures have significant velocity width and spatial elonga- 
tion. Detailed parameters of identified maser features are 
listed in table 2. The X, Y axes and colour show the 
right ascension offset (Aaxcos5), declination offset (AS) 
and Vlsr, respectively. The coordinate origin is (a,<5)j20oo.o 
= (17 h 26 m 01 s .5883, -34° 15' 14". 905), which is the absolute 
position of the brightest maser spot in the second VLBI 
epoch. This position is consistent with CP08 position within 
the range of errors. 

Estimated errors of our position measurements are 
listed in table 1. We have evaluated three error sources, i.e., 
thermal noise in the phase-referenced images, errors of base- 
line ve ctors and atmospheric p hase fluctuations in the same 
way as lNakagawa et aU (|2008h . The atmospheric zenith de- 
lay residuals (~ 0.7 cm on average) are determ ined by the 
image optimizing method (|Honma et al ] 120071 , 2008). The 
largest atmospheric residuals (~ 3 cm) was found in the 
second VLBI epoch, where estimated error reached 1 mas 
in X and 2 mas in Y directions. The relative position un- 
certainty between each maser spot is, on the other hand, 
typically 30 ^tas in X and 100 /xas in Y directions, based on 
thermal noise in the self-calibrated images. 

All maser spots have been found within 200 mas x 100 
mas area. This corresponds to 350 au x 200 au at the source 
distance of 1.7 kpc. We have detected successive maser align- 
ment, which looks like a filament, in the central 50 mas x 50 
mas region. This filamentary structure is pronounced only in 
the flare-up phase tracing shock front (see next subsection). 

Total Vlsr coverage of maser emissions in VERA obser- 
vations is -121 to -5 km s _1 , and almost same as the single- 
dish data. This is well consistent with the blue-shifted part 



of CP08 spectrum. The red end of -5 km s _1 is same as the 
systemic velocity. Extremely red (+87 km s _1 ) component 
reported in CP08 is out of our observing band. Large por- 
tion of flux is concentrated in the range of -80 to -40 km 
s- 1 . 

Spatially, there is no clear gradient of Vlsr for features 
in this dominant velocity range. The features out of this 
range roughly split into east (blue) and west (red) part of 
G353. The most blue and red components are separated only 
100 mas (200 au) along SE - NW direction. The wide range 
of velocity within a very compact region readily dismisses 
the possibility of disc maser which has been propo sed in a 
few cases fe.g jlmai etaL 2006; Trinidad et alj|2009l) . It also 
indicates relatively pole-on geometry of the disc-jet system 
as suggested in CP08. Expected young age of the driving 
source may also be related to such a concentrated distri- 
bution, however, observed proper motions (see next section) 
comparable to the Vlsr values still support moderately pole- 
on geometry even in this case. 



3.2 Flare activity 

We have detected strong maser flare in the second VLBI 
epoch. Time variation of single-dish spectrum after the flare 
is displayed in figure 2. There are two velocity components 
which show significant brightening. One is the component of 
-53 ± 7 km s" 1 (hereafter C50) and another is the compo- 
nent of -73 ± 7 km s _1 (hereafter C70). Almost all emissions 
disappeared in the quiet phase. Figure 3 presents the light 
curves for both flared components during 950 days from Jan 
1st, 2008. The vertical axis indicates peak flux density for 
single-dish data, and sum of the flux densities of maser spots 
detected in the peak velocity channel for VLBI data. Down- 
ward arrows show 3-cr upper limit for C70. The plotted data 
correspond to column 5 of table 3. 

Only C50 had been flared at first and C70 started to in- 
crease its flux about 100 days after observed C50 peak. Sign 
of second flare is detected only for C50 after the 8th VLBI 
epoch (day 841). The systematic increase of the 3-<r upper 
limits in the latter part of C70 is not a signature of sec- 
ond flare, but because of bad weather condition in summer 
season. A variation pattern of H2O maser is, in general, com- 
bination of long-term, globa l pattern and short-term fluctu- 
ation (e.g.. IFelli et al.l 20071 ) . Smoothly varying light curve 
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Figure 1. Spatial and velocity distribution of all maser spots in G353. Colour triangles indicate detected maser spots and their Vlsr 
values. Coordinate origin is 

(«!<5)j2000.0 = (17 h 26 m 0F.5883, -34°15'14".905) which is the position of the brightest maser spot in second 

VLBI epoch. 



indicates that global variation pattern is dominant in C50 
rather than random fluctuation during our monitoring. Al- 
though large dip was caught in day 473 for C70, same global 
pattern was identifiable, i.e., the flare-up continued for ~ 150 
days and rapid decaying. The dip in C70 is roughly factor 
of three weaker than the peak flux density and explainable 
by short-term fluctuation. 

Figure 4 presents time variation of the prominent 
part of the maser alignment that contains most flux of 
C50. This successive alignment had been appeared af- 
ter the first flare. The maximum spatial scale of the 
alignment is roughly 50 mas (~ 90 au) that is al- 
most typical value fo r an individual maser filament in 
an outflow shock (e.g.. | Elitzur. Hollenbach fc McKeelll989l ; 
iMotogi et ai"1l2008l ; iTorrelles et al.ll201lh . A number written 
in upper right part of each panel shows relative day from 
the 1st VLBI epoch (day 183). The alignment was initially 
vertical and dramatically changed to be horizontal during 
the flare phase (2nd - 4th epoch). All maser features in this 
region completely disappeared in the 7th VLBI epoch. An- 
other maser alignment appeared in the final VLBI epoch, 
tracing the second flare. 

Because of non- negligible accelerations of Vlsr, the 
identification of individual maser feature in these alignments 
is difficult (see section 3.3). We, thus, estimated averaged 
proper motions of these alignment, where we averaged the 
positions of all maser spots in each alignment and calculated 
fix and fiy for vertical and horizontal alignments, respec- 
tively. Here, fix and fiy show the linear proper motion in 
each directions. Such a motion perpendicular to a maser 



alignment should trace bulk motion of shock front, if we 
consider the standard model that H2O masers are excited 
in a thin sheet of dense gas compressed by outflow shock 
l|Elitzur et al.|[l989l ). Estimated fix and fiy are 7.8 mas yr _1 
from east to west and 7.6 mas yr~ from north to south, re- 
spectively. These averaged proper motions correspond to ~ 
60 km s _1 at 1.7 kpc and comparable with their Vlsr values 
of -50 - -70 kms" 1 . 

We note that these averaged motions are absolute mo- 
tions respect to the referenced calibrator, and hence, include 
both of the solar motion relative to the LSR and the sys- 
temic motion of G353. The solar motion has already been 
subtracted from averaged motions above. We used the so- 
lar motion estimated from Hipparcos dat a (0.3 mas yr~ for 
X an d -1.1 mas yr" 1 for Y direction. ; iDehnen fc Binnevl 
1998). The systemic motion of G353 is expected to be neg- 
ligible since rotating motion of G353 should be almost par- 
allel to the sun, if we assume that G353 is moving along 
flat Galactic rotation. But if G353 has significant deviation 
from flat rota tion, possible con t ribution can be u p to 15 - 30 
km s _1 (e.g.. iReid et~aI1l2009l ; feaba et al.ll2009r i. However, 
actual systemic motion seems to be enough small because 
no constant and systemic velocity vector has been seen for 
maser features in G353. 



3.3 Maser acceleration 

Variations of Vlsr during our monitoring is plotted in fig- 
ure 5. Plotted Vlsr values are flux density weighted means 
of relevant velocity channels in C50 and C70 (see table 3). 
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Table 2. Detected Mascr Features 



Epoch 


^LSR 


Offset 11 


(mas) 


Total Width 6 


Integrated Flux 




(km s- 1 ) 


X 


Y 


(km s- 1 ) 


(Jy km s _1 ) 


1 


-43.57 


9.65 (0.01) 


42.34 (0.02) 


1.69 


13.01 (0.44) 




-44.21 


10.14 (0.02) 


42.25 (0.04) 


0.42 


0.61 (0.05) 




-46.52 


-0.31 (0.03) 


-0.90 (0.07) 


1.05 


0.89 (0.09) 




-57.9 


-1.65 (0.01) 


15.20 (0.01) 


1.48 


6.19 (0.14) 




-63.38 


9.05 (0.00) 


21.69 (0.01) 


1.90 


7.25 (0.13) 




-66.75 


34.43 (0.01) 


39.46 (0.02) 


1.48 


3.85 (0.13) 




-88.03 


193.03 (0.01) 


66.06 (0.01) 


1.26 


2.70 (0.05) 


2 


-42.69 


1 ^ cn Cn 09l 


k ok Cn (14 s ! 


1.05 


^44 CO 1 (\) 




-43.95 


1 79 CO 01 1 


45 1 7 CO 02) 

t:U. i I i u . Ui y 


1.48 


5 08 CO 1 2) 

U.UO lU.J-il 




-48.8 


47 CO 01 1 

U.rt 1 i u . u _1_ ( 


-ft 95 CO 021 


3.37 


1 34 QQ C1 891 




-51.96 


4 64 CO 001 

l -±.KJ L ± 1U.UU/ 


-5 51 CO 01 1 


5.48 


445 00 C4 09") 




-49.85 


ft 98 (ft 09 s ) 


1 1 4 (ft 03 1 ) 


0.63 


3 78 (ft 1 9 s ) 




-51.54 


ft 1 (S I'D 01 ^ 

U. J-U \\J.\J-L) 


38 (Ci 02 s ) 


2.53 


87 1 9 (^ 28 s ) 

U 1 ■ li 1-L.iOI 




-52.38 




3 14 (ft 01 ") 


0.63 


Q9 (ft ft^W 
u. y i \yj' uu j 




-51.96 


c; 88 fO 00 s ) 

U.OO IU.UU/ 


-7 QQ ("0 01 ") 


0.63 


ft QQ (ft ftft) 

U.£7i? 1 U.UU 1 




-51.96 


-1.10 (0.03) 


3.60 (0.03) 


0.21 


0.10 (0.02) 




-52.17 


-0.07 (0.02) 


1.32 (0.04) 


0.21 


0.61 (0.03) 




-52.38 


-1.12 (0.04) 


4.38 (0.05) 


0.21 


0.07 (0.01) 




-52.59 


-0.22 (0.02) 


2.10 (0.05) 


0.63 


1.41 (0.06) 




-53.43 


-0.55 (0.03) 


3.95 (0.06) 


2.32 


6.29 (0.26) 




-56.17 


-0.79 (0.05) 


7.79 (0.12) 


0.21 


0.01 (0.00) 




-56.38 


-1.19 (0.02) 


14.43 (0.08) 


1.05 


0.29 (0.01) 




-58.7 


-0.39 (0.01) 


18.34 (0.02) 


1.9 


2.35 (0.04) 




-63.76 


107.57 (0.06) 


51.48 (0.11) 


0.42 


0.03 (0.00) 




-66.29 


110.83 (0.04) 


40.11 (0.08) 


0.63 


0.04 (0.00) 




-67.97 


32.99 (0.01) 


43.13 (0.02) 


1.26 


0.63 (0.01) 




-69.87 


33.76 (0.01) 


43.07 (0.02) 


1.48 


0.81 (0.02) 




-72.18 


33.52 (0.02) 


43.85 (0.04) 


0.84 


0.11 (0.01) 


3 


-5.72 


-32 Q3 CO 021 


1 02 08 CO 041 


0.84 


Q4 CO 07) 

U. £7rt I U .UI ; 




-23.63 


-1 1 92 CO 01 "1 

lit \ U . U J- / 


37 09 CO 02) 


1.05 


1 20 CO 051 




-32.68 


-7 70 CO 051 


38 91 CO 061 

tjo.iyj- v u.uuy 


0.63 


1 5 CO 021 




-41.53 


1 9 Q4 fO O 1 ^ 


a fi7 /n 1 3 s ) 


0.21 


09 (ft ftft\ 

U. UZr l^U. UU J 




-43.43 


Q SO TO 04~) 


1-.<-)iJ 1 U.U J J 


1.05 


ft 2Q (ft 01) 




-44.06 


Q 4fi ('O 01 ^ 


44 97 fO 01 ■) 

^±1 .if ^ U . U -L J 


1.48 


9 04 f0 04 s ) 




-45.32 


(S 37 fO 01 ^ 

U. <J 1 \\J . U-L / 


42 1 fi (Y) 03 s ) 

Ti.lU I U.Ud 1 


1.69 


ft ^Q (ft 02 s ) 




-48.06 


3 7fi (ft 01 ^ 


9Q ^8 ("0 01 


2.32 


9 0^ ^O 04 s ) 




-48.9 


-2 00 (0 OS 1 ) 

i . uu \u . Ut_M 


-1 1 ^ TO Ofi'l 

l.li) l U.UU J 


1.26 


ft 38 fO 02 s ) 




-51.85 


9 03 TO 01 'l 


87 fO 03 s ) 


2.95 


37 (ft 9f^ 

U.O / ^U. iU J 




-50.38 


1 2Q (0 01 ^ 


^ ("0 09 s ) 


1.26 


1 27 fO 03~) 

1.^1 IU . UtJ / 




-53.12 


3 39 CO 001 

j.jj \u . uu / 


-7 28 Co no) 

1 .iO 1 u.uu I 


3.58 


1 37 82 CO 401 

lO 1 .Oi lU.rtUJ 




-51.64 


8 74 CO 021 

O. 1 t: I U . Ui ) 


-9 32 CO 02) 

C.ijii lU.UiJ 


2.53 


3 79 CO 1 61 

■J . J C» 1 U. J-U J 




-53.33 


2.71 (0.01) 


-4.75 (0.00) 


0.63 


0.05 (0.01) 




-54.59 


6.44 (0.01) 


-8.69 (0.02) 


1.9 


1.82 (0.05) 




-56.7 


7.06 (0.02) 


-8.31 (0.02) 


1.47 


0.28 (0.01) 




-55.22 


-1.41 (0.03) 


4.03 (0.05) 


0.84 


0.08 (0.01) 




-55.43 


-3.14 (0.04) 


8.13 (0.10) 


0.21 


0.01 (0.00) 




-58.38 


-1.87 (0.01) 


15.95 (0.02) 


1.48 


0.17 (0.01) 




-64.49 


107.28 (0.01) 


50.10 (0.01) 


2.11 


0.39 (0.01) 




-66.81 


108.18 (0.02) 


37.78 (0.04) 


0.84 


0.03 (0.00) 




-66.81 


67.69 (0.04) 


13.31 (0.09) 


0.42 


0.00 (0.00) 




-67.44 


110.58 (0.03) 


38.34 (0.06) 


0.42 


0.00 (0.00) 




-70.39 


32.54 (0.00) 


42.4 (0.00) 


2.74 


1.68 (0.01) 




-72.71 


32.79 (0.00) 


42.81 (0.00) 


2.11 


1.40 (0.01) 




-73.76 


33.00 (0.00) 


42.92 (0.01) 


1.90 


0.58 (0.01) 


4 


-23.91 


-12.31 (0.01) 


34.98 (0.02) 


1.47 


6.83 (0.18) 




-43.93 


83.35 (0.03) 


59.54 (0.06) 


0.21 


0.09 (0.01) 




-44.35 


9.10 (0.06) 


42.72 (0.11) 


0.21 


0.03 (0.01) 




-44.35 


-0.04 (0.01) 


-8.35 (0.02) 


0.84 


1.58 (0.04) 
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Table 2 - continued 



Epoch y LSR 


Offset" 


(mas) 


Total Width 6 


Integrated Flux 


(km s- 1 ) 


X 


Y 


(km s- 1 ) 


(Jy km s _1 ) 


-44.56 


10 CO 01 ") 


-8 42 ("0 02) 


0.63 


73 CO 02*1 


-48.56 


9 (C\ 01 ^ 


^ ^7 (f\ f)°A 

-O'O t \\J.\jO j 


1.48 


9 7fi (f) f\R\ 


-48.77 


^ Sfi (C\ 09 s ) 


97 (0 C\d) 


1.47 


81 (f) 0^^ 


-49.61 


o m (c\ ni\ 

-U.U-L I^U.UO j 


% 01 (f\ C)^} 

-O.UJ- \U.VJO ) 


0.84 


77 (f) 0^ 

U.I/ ^u. UO J 


-52.77 


10 0^ (0 01 1 


-1 9 (0 09") 


1.90 




-51.72 


13.01 (0.02) 


-13.39 (0.05) 


1.05 


0.34 (0.02) 


-53.83 


2.62 (0.01) 


-9.81 (0.03) 


1.69 


2.06 (0.08) 


-53.2 


9.93 (0.01) 


-12.67 (0.01) 


1.05 


1.36 (0.03) 


-55.09 


5.47 (0.01) 


-11.64 (0.02) 


1.05 


1.32 (0.04) 


-54.88 


7.36 (0.03) 


-12.12 (0.03) 


1.26 


0.94 (0.06) 


-56.14 


7.10 (0.02) 


-11.70 (0.02) 


1.05 


0.35 (0.02) 


-56.14 


8 04 (C\ O 1 ^ 

O . \\J . UO J 


1 1 fi1 (0 01} 


1.69 


4fi (0 


-65.62 


107.80 (0.02) 


47.85 (0.05) 


0.63 


0.05 (0.00) 


-69.84 


33.25 (0.02) 


40.05 (0.05) 


0.63 


0.05 (0.00) 


-71.94 


32.86 (0.00) 


41.11 (0.01) 


2.95 


3.10 (0.04) 


-74.47 


33.31 (0.00) 


41.60 (0.00) 


2.74 


14.07 (0.03) 


-95.96 


96.44 (0.01) 


12.48 (0.03) 


0.84 


0.08 (0.00) 


-120.18 


27.84 (0.01) 


8.38 (0.02) 


1.69 


0.39 (0.01) 


5 -43.74 


83.18 (0.03) 


58.33 (0.08) 


0.42 


1.01 (0.13) 


-49.43 


-4.83 (0.01) 


-5.81 (0.02) 


2.11 


11.75 (0.34) 


-54.7 


12.57 (0.02) 


-15.34 (0.06) 


0.42 


0.24 (0.03) 


-56.6 


9.12 (0.03) 


-14.84 (0.05) 


1.05 


0.77 (0.07) 


-75.14 


33 00 TO 01 "I 


41 03 (0 02) 


1.69 


4 67 ("0 1 6"! 


6 -24.13 


-12.84 (0.04) 


34.28 (0.14) 


0.21 


0.25 (0.05) 


-42.88 


-14.26 (0.01) 


1.07 (0.06) 


0.84 


1.36 (0.09) 


-43.72 


83.85 (0.01) 


58.7 (0.03) 


0.63 


0.41 (0.03) 


-45.41 


-5 21 fO 04"l 


-5 68 (0 08") 


0.21 


04 (0 00) 


-47.94 


10.58 (0.06) 


-25.2 (0.09) 


0.21 


0.05 (0.01) 


-49.41 


-4.78 (0.00) 


-5.68 (0.01) 


2.53 


7.10 (0.09) 


-48.78 


-6.67 (0.01) 


1.54 (0.04) 


0.42 


0.08 (0.01) 


-57.42 


9.65 (0.01) 


-14.76 (0.02) 


1.26 


0.54 (0.02) 


-67.32 


72.03 (0.01) 


68.13 (0.02) 


0.84 


0.23 (0.01) 


-75.96 


33.86 (0.00) 


41.62 (0.00) 


1.90 


4.07 (0.04) 


-107.35 


92.64 (0.01) 


38.53 (0.02) 


1.05 


0.39 (0.01) 


7 -50.31 


-5.23 (0.05) 


23.23 (0.11) 


0.63 


0.94 (0.15) 


-56.63 


52.02 (0.03) 


41.85 (0.06) 


0.63 


0.74 (0.07) 


-67.16 


71.42 (0.01) 


69.19 (0.02) 


1.05 


3.96 (0.09) 


-77.27 


33.30 (0.03) 


42.95 (0.07) 


0.63 


0.22 (0.02) 


8 -38.38 


11.95 (0.02) 


-18.89 (0.06) 


0.21 


0.96 (0.09) 


-47.22 


-6.98 (0.01) 


-2.75 (0.03) 


1.26 


6.14 (0.26) 


-47.43 


-25.62 (0.02) 


12.09 (0.13) 


0.21 


0.14 (0.03) 


-47.43 


-20.70 (0.03) 


11.09 (0.13) 


0.21 


0.11 (0.03) 


-48.49 


-5.75 (0.00) 


-4.57 (0.01) 


2.11 


24.19 (0.30) 


-51.65 


-2.10 (0.01) 


-3.84 (0.03) 


1.26 


1.15 (0.05) 


-55.86 


20.96 (0.04) 


-16.42 (0.08) 


0.42 


0.21 (0.02) 


-58.39 


5.95 (0.03) 


-12.31 (0.10) 


0.21 


0.03 (0.00) 


-67.02 


73.96 (0.02) 


69.02 (0.06) 


0.42 


0.08 (0.01) 



a Relative offsets from (a,<5)j 2 ooo.o = (17 h 26 m 01 s .5883, -34°15'14".905). 
b Full width at zero intensity (FWZI) for each maser feature. 



All parenthetic values indicate associated errors. 



Associated error bars show standard deviations from esti- 
mated means. These are good indicator of how maser emis- 
sions concentrate in the velocity domain. Filled squares and 
triangles correspond to single-dish and VLBI data, respec- 
tively. We also distinguish different peak flux densities by 
colour coding. Linear acceleration is evident in C70 during 
all detected epochs. 



Systematic acceleration is also seen for C50 only in the 
flare-up phase (day 300 - 600), although these data points 
have relatively larger deviations compared to C70. Linearly 
fitted acceleration is -5.4 ± 1.2 km s _1 yr _1 . This is con- 
sistent with the value of -4.8 ± 0.5 km s _1 yr for C70 
which is isolated from other components in both of velocity 
and spatial domains. The fact that two independent compo- 
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-70 -60 

%R [ km s_1 ] 



Figure 2. Time variation of singlc-dish spectra. The vertical axis is in a hundred Jy unit. The numbers on far right show the relative 
day from Jan 1st, 2008. The systemic velocity of -5 km s _1 and several weak high velocity components (< -100 km s _1 ) arc out of the 
plot range. 



C50 Light Curve 




150 



100 



C70 Light Curve 
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450 600 
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450 600 
Day 
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Figure 3. The light curve for C50 (left) and C70 (right). Filled squares and triangles correspond to single-dish and VLBI data, 
respectively. The horizontal and vertical axis shows the relative day from Jan 1st, 2008 and the peak flux density, respectively. The latter 
is sum of the flux densities of all detected spots in the peak velocity channel for VLBI data. Downward arrows present 3-cr upper limit. 
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Figure 4. Time variation of the prominent maser alignment around the brightest maser spot. Colour 
triangles indicate detected maser spots and their Vlsr- The numbers and associated parenthetic values in 
upper left side of each panel show the day offset from the first VLBI epoch and corresponding observing 
epochs, respectively. Two labelled velocities show estimated average proper motions assuming the distance 
of 1.7 kpc. 
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Day 

Figure 5. Variations of the flux density weighted mean of Vlsr 
for C50 (top) and C70 (bottom). The definitions of horizontal 
axis and markers are same as figure 4. The colour scale shows the 
peak flux densities. Error bars indicate standard deviations from 
the weighted means, which are measure of flux concentrations in 
the velocity domain. 

nents have shown same accelerati on suggests that it is not 
apparent acceleration discussed in lBrand et al.l l)200ot ). but 
real one caused by a momentum suppl y from a protostella r 
outflow or wind as is the case of W75N l|Hunter et al.lll994 ). 

In addition to the acceleration, two discontinuous jumps 
of Vlsr occurred in C50. The first jump was seen at the ini- 
tial flare event (~ day 200) and the other one was synchro- 
nised with the run-out of the horizontal alignment (~ day 
600). The initial jump clearly trace the flare-up of the maser 
alignment and the second jump is caused by attenuation and 
extinction of higher velocity (~ -56 km s _1 ) components in 
C50. 



4 DISCUSSIONS 
4.1 Outflow activity 

iBrand et all (|2003h has attributed one possible cause of 
long-term, grobal variation of H2O masers to a variation of 
outflow activity itself. This seems to be applicable for G353, 



because it is difficult to fully attribute the systematic change 
of the shock front to internal fluctuations of coherent path. 
Synchronised acceleration of C50 and C70, which are sepa- 
rated ~ 100 au in space and 20 km s _1 in velocity, implies 
that subsequent momentum supply accelerates entire part 
of masing gas sheet, keeping or gradually losing coherence. 
If this is the case, intermittent maser flare of G353 is simply 
interpretable as episodic launching of an outflow. Expected 
launching time-scale is only ~ 2 yr based on the interval 
between two flares in this case. This is well comparable wi th 
the typical time-scale of several radio jets (|Angladalll996h . 

The characteristic velocity of ~ 500 km s _1 , which is 
estimated from the spatial scale of maser distribution (~ 
200 au) divided by the time-scale (~ 2yr), is also consistent 
with the proper motion veloci ty of radio jets associated with 
several high mass YSOs (e.g.. lMarti. Rodriguez fc Reipurthl 
1 19981 ; ICuriel et al.1 12006). Thus, observed variation can be 
caused by standard radio jet in high mass star-forming re- 
gion. The relation between the expected velocity of a jet and 
the 3D velocity of masers (~ 10 km s" 1 ) is roughly same as 
the case of CepA HW2 region l|Torrelles et al.ll201ll ). where 
velocities of radio jet and relatively wide angle maser flow 
is 500 and 70 km s" 1 , respectively. G353 is thought to be a 
possible candidate of high mass YSO with a radio jet activ- 
ity through these reasons. 

The momentum rate of maser flow is estimated as 1.1 x 
10 -3 Mq km s _1 yr -1 based on the acceleration of ~ 5 km 
s _1 yr -1 . Here, we assumed that the masing gas sheet has 
uniform thickness of 1 au, i.e., typical size of maser feature. 
Total area of the sheet is evaluated as 200 x 200 au 2 from 
the extent of maser distribution. Adopted number density 
of molecular hydrogen is 10 9 cm -3 that is also typical value 
for a H2O maser feature. We note that assumed volume is, 
and hence, the momentum rate is only a lower limit, be- 
cause a mass fraction of masing gas shee t to total shocked 
gas should be small l)Elitzur et al.l Il989l ) . Estimated lower 
limit is clearly larger than a typical v alues of low mass star - 
formation (~ 10" 5 M & km s" 1 vr -1 ; lArce fc Sargentj|2006l ) 
and consistent w ith that of high mass star-formation (e.g., 
lArce et "aJ]|2007l for comparison). If we adopt a jet velocity of 
500 km s _1 , it is immediately converted to the outflow rate 
of 2.2 x 10 -6 Mq yr -1 . We finally get the lower-limit accre- 
tion rate of 10 -5 - 10 -4 Mq yr -1 , where we have adopted 
typical fraction between outflow rate ove r accretion rate (1 - 
10 % ) in a case of MHD dri ven flow (e.g JPelletier fc PudritS 
ll992l ; lMachida et al.|[2006T ). 

Although we have considered pu re MHD jet (e.g., 
iMachida. Inutsuka fc Matsumotol l2008h here, strong radia- 
tion from the host YSO should also play an important role 
for outflow launching in high mass star-fomation. If the host 
YSO is enough young and an accretion rate is enough high, 
accretion luminosity i s thought to be sign ificant or even 
dominant in total flux (|Hosokawa et al.ll2010l ). Outflow vari- 
ations may directly reflect variations of accretion rate onto 
the protostellar surface in this case. If this is true, varia- 
tion pattern of the H2O maser can be important clue for 
accretion mechanism at the innermost part. Further discus- 
sion about this requires exact stellar parameters in IR wave- 
length. 

Another explanation for intermittent maser flare is a 
precessing jet, where we can observe maser emission only in 
a certain fraction of precession period because of beaming 
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Table 3. Parameters of the maser components 



Day 




AV b 


Total Width Peak Flux Density 


Integrated Flux 




(km s 1 ) 


(km s _1 ) 


(Kin s j 




(Jy km s _1 ) 








C50 single- 


dish data 




426 


-53.20 


1.63 


8.40 


174.67 (3.74) 


R1 7 Q1 ('\ 7A\ 
ui / ■ y ± ii. ill 


434 


-53.49 


1.48 


5.67 


138.75 (8.01) 


4XQ 77 f3 fi7*l 


437 


-53.43 


1.48 


5.67 


122.79 (6.54) 


41 9 O^i f3 00*1 


442 


-54.00 


1.70 


5.88 


62.63 (5.27) 


31 9 77 (1 40 "l 


473 


-54.94 


1.40 


3.99 


64.57 (8.01) 


904 1 Q (1 f\7\ 


500 


-54.13 


1.67 


4.62 


68.64 (6.14) 




501 


-53.86 


1.77 


3.78 


76.39 (7.41) 


91 9 (1 41 ^ 


551 


-55.36 


0.85 


0.42 


40.73 (11.95) 


93 08 C5 48 "l 


671 


-48.86 


0.27 


0.63 


24.30 (5.74) 




686 


-50.97 


3.49 


1.05 


28.38 (4.74) 


4^ i q (o 1 d) 


702 


-48.81 


0.56 


1.89 


14.56 (2.87) 




716 


-48.64 


0.55 


1.68 


12.29 (2.27) 


^9 Q9 (1 07 s ! 


729 


-50.00 


2.96 


1.26 


10.88 (2.47) 


94 ^7 M 1 4^ 


744 


-48.81 


0.74 


0.84 


10.55 (2.80) 


90 ^7 (1 97 1 ) 


759 


-48.71 


0.66 


1.47 


1 5 62 (1 47*1 


Qfi 7Q (1 14^ 


780 


-49.53 


3.24 


1.89 


17.76 (3.27) 


48.41 (1.47) 


796 


-48.44 


0.31 


0.84 


23.44 (5.14) 


66.10 (2.40) 


832 


-49.34 


2.65 


5.67 


96.08 (2.80) 


204.39 (1.34) 


834 


-49.15 


2.30 


6.09 


101.63 (3.27) 


211.07 (1.47) 


839 


-49.35 


2.79 


6.72 


114.65 (3.14) 


250.66 (1.47) 


852 


-48.68 


2.06 


4.41 


108.64 (4.27) 


235.30 (1.94) 


874 


-48.11 


1.37 


2.94 


151.04 (8.15) 


283.25 (3.74) 


893 


-48.27 


1.82 


3.15 


218.41 (9.15) 


368.71 (4.21) 


899 


-47.95 


0.46 


1.89 


125.33 (11.62) 


210.26 (5.34) 


901 


-47.98 


0.48 


1.89 


154.31 (15.96) 


275.03 (7.34) 


906 


-47.82 


0.56 


2.31 


153.64 (11.42) 


272.76 (5.27) 


912 


-47.93 


0.31 


1.05 


142.02 (26.71) 


234.97 (12.29) 


913 


-47.88 


0.43 


1.68 


189.97 (17.09) 


310.15 (7.81) 


C50 VLBI data 


183 


-56.65 


3.57 


2.52 


23.50 (0.99) 


24.00 (0.56) 


323 


-51.24 


1.70 


12.81 


716.30 (17.01) 


2224.24 (14.28) 


402 


-52.93 


1.15 


14.07 


632.37 (2.86) 


969.45 (3.29) 


503 


-53.13 


2.50 


9.24 


39.42 (1.39) 


122.53 (1.51) 


624 


-51.08 


2.83 


4.20 


12.80 (0.94) 


21.95 (0.12) 


674 


-50.30 


3.13 


4.62 


11.26 (0.40) 


23.36 (0.05) 


757 


-54.61 


2.99 


1.26 


2.17 (0.37) 


2.10 (0.07) 


841 


-48.94 


2.19 


6.51 


66.42 (0.63) 


87.55 (0.12) 








C70 single- 


dish data 




426 


-73.79 


1.24 


4.20 


106.50 (3.74) 


203.39 (1.94) 


434 


-74.13 


0.71 


2.94 


122.39 (8.01) 


233.3 (4.14) 


437 


-74.12 


0.75 


3.15 


117.25 (6.54) 


249.53 (3.34) 


442 


-73.88 


0.87 


3.36 


82.73 (5.27) 


202.32 (2.67) 


473 


-74.26 


0.41 


1.47 


49.74 (8.01) 


95.68 (4.14) 


500 


-74.45 


0.34 


1.47 


151.44 (6.14) 


168.47 (3.14) 


501 


-74.31 


0.72 


1.89 


155.98 (7.41) 


195.24 (3.81) 


551 








< 35.86 (11.95) 




671 


-75.95 


0.20 


0.63 


23.44 (5.74) 


35.12 (2.94) 


686 


-76.28 


0.25 


0.84 


41.87 (4.74) 


55.75 (2.40) 


702 


-76.72 


0.24 


0.84 


19.16 (2.87) 


11.02 (1.47) 


716 


-76.72 


0.25 


0.84 


9.68 (2.27) 


3.54 (1.14) 


729 


-76.81 


0.11 


0.21 


7.34 (2.47) 


7.41 (1.27) 


744 








< 8.41 (2.80) 




759 








< 7.41 (2.47) 




780 








< 9.82 (3.27) 




796 








< 15.42 (5.14) 




832 








< 8.41 (2.80) 




834 








< 9.82 (3.27) 




839 








< 9.41 (3.14) 




852 








< 12.82 (4.27) 
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Table 3 - continued 



Day Vlsr" A V b Total Width Peak Flux Density Integrated Flux 

(km s -1 ) (km s" 1 ) (km s _1 ) (Jy) (Jy km s — 1 ) 

C70 single-dish data 



874 
893 
899 
901 
906 
912 
913 



< 24.44 (8.15) 

< 27.44 (9.15) 

< 34.85 (11.62) 

< 47.88 (15.96) 

< 34.25 (11.42) 

< 80.13 (26.71) 

< 51.28 (17.09) 



C70 VLBI data 



323 


-69.44 


1.40 


4.41 


16.60 (0.68) 


38.03 (0.51) 


402 


-71.92 


1.67 


7.35 


18.90 (0.22) 


65.93 (0.28) 


503 


-74.04 


1.14 


6.51 


165.00 (0.26) 


175.37 (0.44) 


624 


-75.41 


0.38 


1.89 


18.60 (0.80) 


19.77 (0.23) 


674 


-75.99 


0.38 


2.10 


28.80 (0.15) 


24.79 (0.23) 


757 


-77.32 


0.25 


0.63 


2.04 (0.34) 


1.08 (0.11) 



a Flux density weighted mean. 

b Standard deviations from the weighted means. 

All parenthetic values indicate associated errors. 



effect. If simple binary interaction which has rotating period 
of 2 yr is considered, expected binary separation is over 3 
au with a total mass of 10 Mq. By contrast, if we assume 
escape velocity at he launching point as typical velocity of 
a jet, expected launching radius is estimated as follows, 

here M bj ec t and Vjet indicate mass of the driving source and 
velocity of launched jet, respectively. It is well smaller than 
the expected binary separation and this scenario seems to be 
possible. High-resolution image of molecular lobe obtained 
with upcoming Atacama Large Millimeter / submillimeter 
Array (ALMA) will be helpful to distinguish two scenar- 
ios. Direct detection of precessing trajectory of a jet is also 
the best way to test them. This can be achieved by high- 
sensitivity observation of radio continu um or IR line emis- 
sions, as the case of IRAS 16547-4247 l|Brooks et al.l 120031 ; 
iRodrfguez et al.ll2005l ). 



4.2 Driving source 

The H2O maser distribution is completely inside the error 
circle of the position of class II CH3OH maser reported by 
CP08. Since there is no d etectable OH maser emission (e.g., 
CP08; iBreen et alj|20ich . these masers are thought to be 
excited by identical high mass protostellar object. If the H2O 
masers in G353 is actually excited by a radio jet, such a close 
location (< a few hundred au) of two types of masers and 
the radio jet is also similar to the cases o f several high mass 
proto stellar objects. Besides CepA HW2 (|Torstensson et all 
l20Tll ). IR AS 20126+4104 is the best example of th is kind 
of source (Hofn er et al. I l2007l ; iMoscadelli et al.|[201ll ). 

H 2 masers in IRAS 20126+4104 are obviously ex- 
cited in front of the rad io knot which emerges from the 
nearly edge-on disc (e.g., lEdris et al.ll2005l ; ICesaroni et al.l 



120051) . The intriguing feature of IRAS 20126+ 4104 is asym- 
metry of the radio knot and H2O maser in I Hofner et al.l 
(|2007l ). In spite of overall bipolar structure of outflow ver- 
ified by radio and IR line observations dSu et alj 120071 : 
ICaratti o Garatti et al.ll2008l ). the radio knot and the most 
of H2O masers have been detected only in blue-shifted side. 
The non-negligible blue-shift dominance of H2O maser is 
also evident in long-term single-dish study (|Felli et al.ll2007| ) . 
Their light curve, furthermore, has shown weak intermittent 
variation of 2 - 3 yr scale alike G353. 

These facts strongly support the interpretation that 
G353 is a candidate of high mass protostellar object with 
intermittent maser activity caused by a protostellar jet, as 
proposed in CP08. If this is the case, G353 can be a good 
target of direct imaging of the accretion disc, because of 
the possible pole- on geometry of the disc-jet system (e.g., 
iKraus et alj|20ich . 



5 SUMMARY AND CONCLUSION 

Our VLBI and single-dish monitoring of 22 GHz H2O 
maser emission from the high mass young stellar object 
G353. 273+0.641 with VERA and Tomakamai 11-m radio 
telescope has detected two times of maser flares. The eight 
epochs of VLBI data have revealed that these flares have 
been accompanied by structural change of the prominent 
shock front traced by H2O maser alignments. We have de- 
tected only blue-shifted maser emissions and all maser fea- 
tures have been distributed within very small area (200 x 
200 au 2 ) in spite of wide velocity range (> 100 km s _1 ). 
These facts suggest relatively pole-on geometry of the disc- 
jet system. 

The light curve shows notably intermittent varia- 
tion. This suggests a possibility that the H2O masers in 
G353. 273+0. 641 are excited by episodic radio jet. The inter- 
val of two flares is about two years and well comparable to 
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typical time-scale of a radio jet. The characteristic velocity 
of 500 km s -1 , which is estimated from the time and spatial 
scale, is also consistent with the velocity of the radio jets as- 
sociated with several high mass young stellar objects. These 
facts also suggest that the H2O maser in G353. 273+0. 641 
is excited in a strong shock of a radio jet. The precessing 
jet may be also able to explain intermittent maser flare and 
this will be tested by high-resolution imaging of molecular 
lobe or direct detection of precessing trajectory of a jet. 

Two strong velocity components of C50 (-53 ± 7 km 
s _1 ) and C70 (-73 ± 7 km s _1 ) have shown synchronised 
linear acceleration of the flux weighted Vlsr values (~ -5 
km s~ yr~ ) during the flare phase. The acceleration can 
be converted to the lower-limit momentum rate of 1.1 x 
10~ 3 M km s" 1 yr _1 . If we use a jet velocity of 500 km 
s _1 , expected lower limit of accretion rate is 10 -5 - 10~ 4 
Mq yr^ 1 in the case of MHD driven jet. 

The co-existence of class II CH3OH maser at 6.7 GHz 
and absence of any detectable OH masers, as pointed-out in 
CP08, indicate that the driving source is a high mass proto- 
stellar object. If H2O masers are actually excited by a radio 
jet, the close location (< a few hundred au) of two types of 
masers and the radio jet is also similar to the cases of several 
high mass protostellar objects. IRAS 20126+4104 especially 
has H2O masers excited by the radio jet and quite resembles 
G353. 273+0. 641 in blue-shift dominant spectra and short- 
term (2-3 yr) intermittent variation of the maser. We con- 
clude that G353 is a candidate of high mass protostellar 
object alike IRAS 20126+4104, as proposed in CP08. The 
possible pole-on geometry of disc-jet system can be suitable 
for direct imaging of the accretion disc in this case. 
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